COMPUTING SL(2,C) CENTRAL FUNCTIONS WITH 

SPIN NETWORKS 



SEAN LAWTON AND ELISHA PETERSON 

Abstract. Let G ~ SL(2,C) and be a rank r free group. 
Given an admissible weight A in N"^''^'^, there exists a class func- 
tion defined on Hom(Fr,G') called a central function. We show 
that these functions admit a combinatorial description in terms 
of graphs called trace diagrams. We then describe two algorithms 
(implemented in Mathematica) to compute these functions. 



1. Introduction 

Let Fj. = (xi,...,Xr) be a rank r free group and let G = SL(2,C). 
The representations 9^^ = Hom(Fr, G) are an affine variety isomorphic 
to G^''. G acts on 9^ by conjugation. 

As an affine variety is associated to a reduced C-algebra C[lHr], 
and as G-modules we have the following decomposition: 

|fc|=lJl 

where Vk = Sym'^(C^), z = (ii, Z2, ir) G , \i\ = h + ■ ■ ■ + ir, and 
j = (ji, ...,jr-i) G N*"^"*^ satisfies j G \i\ if and only if 

0<ji< min(zi H \- ii - 2{ji H h ji-i), ii+i) 

for all / = 1,2, ...,r- 1. 

This generalizes the complexification of the Peter- Weyl Theorem 

feeN 

Definition 1.1. Given the above isomorphism, for each triple i,j,k 
such that i G N^', j, k E \i\, and \j\ = \k\, there exists a class function 
X'{'" G End(V^|^|_2|j|))'^ C CfD^r]*^- We refer to the functions Xi'* as rank 
r central functions. 
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Given any two central functions x ^i-nd x', their product xx' is again a 
polynomial function in C[lKr]'^, and thus xx' = Yl ^kXk for coefficients 
Afc and other central functions Xk- This ring structure is very far from 
the natural ring structure of C[!Hr]'^ C C[9^r], and is quite mysterious. 
It seems, thus far, to be best understood in terms of special graphs 
called trace diagrams. 

In |LP09j we described rank 1 and rank 2 central functions and deter- 
mined their ring structure by exploring the calculus of trace diagrams. 
It is the purpose of this paper to show that this association generalizes 
and to describe algorithms to compute these functions explicitly in the 
case of rank 3 central functions. 

The diagramatic form, shown in Figure [H can be deduced to exist 
from work of Baez |Bae96] and later Sikora [SikOlj . One of the most 
important points of this paper is that central functions can be under- 
stood entirely in terms of combinatorics given by a graphical calculus. 
This culminates with Theorem 14. 10[ 



Figure 1 . Diagrammatic Form of Central Functions 

The paper is organized as follows. In the next section, we introduce 
spin networks and trace diagrams. Then we go through the construc- 
tion of rank r central functions and work out a few examples in section 
3. In section 4 we further the study of the trace diagram calculus 
in preparation for the algorithm to compute them. In section 5 we 
give two algorithms to compute rank 3 central functions that we have 
successfully implemented in Mathematica. In the appendix, we review 
some relevant invariant theory and representation theory. 

2. Spin Networks and Trace Diagrams 

The tools used in this paper to explore central functions are trace 
diagrams, which are a slight generalization of spin networks, which will 
be formally defined in Definitions 12.11 and 12.21 

Informally, a spin network is a graph that can also be interpreted 
as a function. Depending upon the type of diagram, the function may 
be between tensor products of V = C^, or between tensor powers of 
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SL(2, C) representations, which are symmetric powers of V. In order 
for a spin network's function to be well-defined, the inputs and outputs 
of the function must be specified, and the graph must be given a small 
amount of extra structure. Once this is done, algebraic rules become 
local diagrammatic rules that are much easier to work with. 

Trace diagrams generalize spin networks by allowing matrices to act 
on specific strands of the diagram. When the diagram has no inputs or 
outputs, it becomes a function SL(2, C) x ■ ■ ■ x SL(2, C) — C that is 
invariant under simultaneous conjugation, hence an element of C[lHr]'^ 
for some r. 




Figure 2. Example of a spin network representing a 
function from Vi^Va to V2® V3®V4, and a trace diagram 
representing a function SL(2, C) x SL(2, C) — > C. 

Diagrammatic notations such as spin networks were developed by 
physicists in the middle of the 20th century as a graphical descrip- 
tion of quantized angular momentum. The theory has been developed 
extensively since that time, most notably by Cvitanovic |Cvi08j and 
Stedman [SteQOj . Trace diagram calculations are very similar to skein 
module calculations in knot theory. Indeed, the Kauffman bracket 
skein module can be thought of as a quantization of the diagrams in 
this paper |Kau91] . 

This section follows the diagrammatic conventions of |LP09t IPetOGt 
IPet09j . rather than the more common conventions of [CFS95t IKau91t 
Maj99] . The most obvious difference is (12.41) , which means that several 
diagrams will differ in sign from these sources. Proofs are omitted, as 
they do not contribute to the main point of the paper. We refer the 
reader to |LP09t IPetOGl IPet09j for these. Similar proofs for alternate 
diagrammatic conventions are in |CFS95[ IKau91] . 

2.1. Diagrammatic Description of V^. We represent vectors v G 
y = by |j , and dual vectors v* E V* hy ] . Strands without 
vector labels are unknown vectors. Tensor products are represented 
by drawing diagrams adjacent to each other: f 1 (S) f 2 = J-J-,. The 
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symmetrization operation is indicated by 

JJ L±, 1 X ^ JJ L±, 

(2-1) ^=;^E^' 

which is an idempotent since 

JJ Li, JJ L±, JJ L±, 

I n I o I n I = I n \ 

II ■■■ II II--- II II-- - II 

The basis elements for Vn = Sym"(y) and V* = Sjm^{V*) described 
in section lA.ll take the form 

n — k k 



1 n 1 


1 

- - 


1 1 

- [£ij l£2j - - 


1 

- l£2J 



------ - l' n ' ^ 

1 r 

It is also customary to represent elements of Vn and V* by thicker 
strands labeled by n, so that the basis elements just described may also 
be written 

(2.3) 



H - 

[tr 



The pairing between Vn and V* is given by 

[n^rri 

= Ski I 



0- 



The SL(2, C)-action on Vn is represented by a matrix g G SL(2,C) 
drawn on the strand: 

■ n 



9 ■ Hn-fe = 



2.2. Spin Networks. The formal definitions of spin networks and 
trace diagrams follow. 

Definition 2.1. A spin network is a trivalent graph drawn in the plane 
with edges labeled by representations o/SL(2, C), such that the labels 
adjoining each vertex form an admissible triple. The graph is drawn 
with all free ends at the top or the bottom, and all local extrema are in 
general position relative to the orientation of the diagram. 

Definition 2.2. A trace diagram is a spin network with one or more 
edges marked by elements of the matrix group SL(2, C), with no marking 
occurring at a local extrema. 
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Every trace diagram corresponds to a unique function. The domain 
of the function is a tensor product of irreducible SL(2, C) -represent at ions, 
given by the labehng of the input edges, while the co-domain is given 
by the labeling of the output edges. 

The function may be computed by decomposing the diagram into 
horizontal slices. The diagram's function is the composition of the 
functions of the individual slices. (One can also decompose the diagram 
into vertical slices, which are conjoined as tensor products.) 

We now describe the lexicon for reading off a spin network's function. 
Let {61,62} be the standard basis for C^. The single-strand cup and 
cap diagrams are 

(2.4) V J = 61 (8> 62 — 62 (8> 61 and : v ® w det[v w], 

respectively. Building upon this definition, the symmetrized cup and 
cap maps are 

(2.5) \^-^[^ and ^ 

From (12.41) and the action of SL(2, C) on = C^, one can deduce 

(2.6) (y) = tT{g). 

This is why the diagrams are referred to as trace diagrams. 

Using Schur's Lemma, the Clebsch-Gordan decomposition ( lA.ll) im- 
plies that for each admissible triple {a,b,c}, as specified in Definition 
1X31 there are unique (up to a nonzero multiple) injections Vc ^ Va^Vf, 
and projections K ® H -» V^. These are sometimes called intertwining 
operators. The injection has the following diagrammatic form: 





(2.7) 

c 'T n r 

The projection (8> H ^ K; is obtained by vertical reflection. 
The following notation will be used extensively in later sections. 

Notation 1. Given an admissible triple {a,b,c}, define 

. b + c - a /,x a + b + c 
^a{b, c) = , e(a, b, c) = . 

In the figure ta{b,c) represents the (unique) number of strands con- 
necting the b and c edges in the expansion. In these terms, the admis- 
sibility condition in Definition IA.3I can be reformulated as 

taib, c) > 0, Zbia, c) > 0, ec(a, b) > 0. 
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2.3. Spin Network Relations. Spin networks satisfy certain skein 
relations that can be leveraged to reason about the underlying func- 
tions. The spin network skein module is the space of formal sums of 
spin networks with coefficients in C, modulo certain relations arising 
from the representation theoretic "building blocks" fl2.ip -( l2T71) . 

The relations that follow can be proven directly from spin network 
definitions. Proofs are contained in |LP09] and earlier works. 

A general spin network skein relation has the form X^m'^m^m ~ 0' 
where ^ C and each s^, represents a diagram. Such relations are 
well-behaved under refiections: 

Proposition 2.3 (Spin Network Refiection, Proposition 3.6 in |Pet06] ). 

Given a spin network s, denote by si the spin network obtained by ver- 
tical reflection ofs, and by s*^ the spin network obtained by horizontal 
reflection ofs. Then 

^ ttmSm = <^ araSl, = <^ ^ ft^S;;; = 0. 

m m m 

Other kinds of topological moves may introduce signs. This is the 
case for local extrema introduced into a strand, crossings of edges ad- 
jacent to a vertex, and reorientation of vertices, as indicated in Figure 

B 



Figure 3. Diagrammatic moves on single strands in a 
spin network that introduce signs. 

Define Sa{b,c) = (-1)^«(M. 

Proposition 2.4 (Spin Network Sign Changes, Proposition 3.22 in 
[LPUQ] ). 



(2.8) = (-l)f 



(2.9) 



X =5,{a,b) A =5„(6,c) 



A more general version of this result follows. See Figure [3] for the 
meaning of 'kink' and 'crossed extrema'. 

Proposition 2.5. Let Ti and T2 be topologically equivalent trace dia- 
grams. Then 

Ti = (-n'^+^T^, 
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where k is the number of "kinks" andx the number of "crossed extrema" 
in the diagrams obtained by expanding edges as in f l2.7l) . 

The following two propositions describe the most basic diagrams, as 
well as how to join and separate strands in spin networks. 

Proposition 2.6 (Closed Spin Networks, Proposition 3.19 in |LP09j ). 



= A(c) = c+l = dim(K); 

Caib, c)\tb{a, c)!ec(a, 6)!(e(a, 6, c) + 1) 



b,c) 



albld 



Proposition 2.7 (Bubble/Fusion Identities, Propositions 3.20,3.21 in 
[LPMl). 



(2.10) "(^^ = be(a,6)|5,,, 

a. jb °W'' 

(2.11) = E X ' 

w;/iere b,(a,6) = = 

We refer to bc(a, 6) and fc(a, 6) as the bubble constant and fusion 
constant, respectively. 

2.4. SL(2, C)-Equivariance and Trace Diagrams. The rules for 
manipulating trace diagrams are the same as those for manipulating 
spin networks, since all parts of a diagram, as mappings, are SL(2, C)- 
equivariant. In particular: 

Proposition 2.8 (Equivariance of Trace Diagrams, Propositions 3.8 
and 3.11 in [LP09j). For g e SL(2, C), 



r'-' '-^ 

n 



n 

TT — — rr 



(ff) (5) 





a) 
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Note the consequences 

rj^U^ and U"^^. 

These equivariance properties permit every relation of the previous 
section to be adapted to include matrices along the edges. In later 
sections, we will often display diagrammatic manipulations without 
the required matrices, appealing to equivariance to show that they also 
hold with the matrices. 



3. Rank r Central Functions 

3.1. Decomposition Theorem. Let C[SL(2,C)] be the coordinate 
ring of the variety G = SL(2,C), and recall our convention that Vk 
denotes the k^^ symmetric tensor of C^, denoted by Sym'^(C^). 

The following theorem is a consequence of the "unitary trick", the 
Peter- Weyl Theorem, and the fact that the set of matrix coefficients of 
SL(2, C) is exactly its coordinate ring. See |LP09] for a detailed proof. 

Theorem 3.1 (Decomposition). There is an SL{2, C) -module isomor- 
phism 

C[SL(2, C)] ^ 5^ ® \4 ^ 5^ End(Vfc). 

keN fceN 

The isomorphism is given by defining 

T:^K*®K^C[SL(2,C)] 

n>0 

by linear extension of the mapping 
where X = ( "^^^ ] is a generic matrix. 



X21 X22^ 

In particular, letting {61,62} be the standard basis for C^, 



n*_;.(X ■ n„_,) = ® = n*_fc ((xii6i + X2162)" '(a;i26i + X2262)') 



(3-1) = E ©"'rroo 



^ (n~l\(l\ n-l~i l-j i j 

jj-^ii X12 X21X22- 



i+j=k 
0<i<n-l 
0<j<l 
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3.2. Applying the decompostion. 

C[SL(2,C)^n ^ C[SL(2,C)]®" 

- (8) 

l<k<r VifeeN / 

(n,...,jr)eN'- 

~ 5Z v;* ® ■ ■ ■ ® v;* (g) (g) ■ ■ ■ (g) . 

(n,...,ir)GN'- 

As is stated in |Pet06] . the isomorphism above is determined by the 
following association: (e*^ (g e*^ (g ■ ■ ■ (g e*^ ) (g (e^j (g (g ■ ■ ■ (g e^v) maps 
to the polynomial function 

(Xi,X2, ...,X,) ^ e;(Xi ■ e,Je*,(X2 ■ e,J ■ ■ -e^X, ■ e.J. 

We will call this "tensorial contraction." 

Our principal interest is with the invariant polynomial functions that 
arise in this fashion. To determine these polynomials we will need a 
notion of "admissibility." 

We say {{Vi-^^, ...,Vi^},Vx) is an admissible pair if and only if \4 
Vii (g ■ ■ ■ (g is a G-module injection. In this case there exists a G- 
module so (g ■ • ■ (g V^^ ~ © (as G-modules). The existence 
of an injection corresponds to a way to connect a single x strand to the 
ii,i2, . . . ,ir strands in an admissible way: 




Several injections are possible, but in this paper we focus on the left- 
associative injection 



(3.2) 




This diagram is only admissible if the triples at each vertex are admis- 
sible, meaning mi G [ii,i2j, and for / > 1, m/ G [mj_i, ii+ij . 

From Proposition IA.2t {V^^, Vi^} and are admissible if and only if 
X = ii + i2 — 2ji for < j i < min(ii, 12). 
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Now consider {V^j, V^j, V^g} and V^. Using the above example and 
Proposition IA.2I a second time we have 



where < ji < min(ii, 22) and < j2 < min(zi+i2— 2ji, 23). Therefore, 
({^1! 5 ^ia}' is an admissible pair whenever x = ii + i2 + is — 
2(ji + J2) and both inequalities < ji < min(?i,Z2) and < ^2 < 
min(zi + 12 — 2ji, ^3) are satisfied. 

Generalizing these examples by iteratively using the Clebsch-Gordan 
formula to decompose V^^ ® ■ ■ ■ we come to the following notation 
and definition. 

Let i = {ii, 12, ir) £ N*", and let \i\ = ii + ■ ■ ■ + ir- 

Definition 3.2. We say that j = {ji, jr-i) G W^^ is i- admissible 
{and denote it by j G \i\) if and only if for all \ <l <r — 1 we have 

0<ji< min(zi H \- ii - 2(ji H h ii+i). 

Note that this is precisely the condition mj G \ii,i2\,fni & \^i-i,ii+i\ 
given earlier with mi = ii -\ h ii+i - 2(ji H \- jn)- 

Recall G = SL(2,C). We then use Clebsch- Gordon iteratively with 
respect to Theorem 13.11 to conclude: 



'i+i2- 



2n ®Vi,^Y^ Vi 



:i+«2+i3-2(il+i2) 



J2 V*^---(E)V*(g)Vi,(^-- 



(n,...,ir)GN'- 





Since the above maps are SL(2, C)-equivariant, 



) 



G 



2\k 



By Schur's Lemma, 




Therefore, 
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2|j| 
\k\ = \l\ 

Definition 3.3. Given the above isomorphism, for each triple i,j,k 
such that i , j,k ^ \i\, and \j\ = \k\, there exists a class function 
Xi''' G C[9^r]*^ which corresponds to a generating homothety {unique up 
to scalar) in End(V^|j|_2|j|p'~^- refer to the functions Xi''° as rank r 
central functions. 

Denote by CXp C C[d\rf the linear span over C of Xp. 
In these terms, 

\k\ = \j\ 

Thus, the central functions Xi''' form an additive basis for C[9^r]'^- 
However, the multiplicative structure in terms of this basis is very 
complicated and not at all obvious. 

We note that i has r entries, k and j have r — 1. So the index relation 
\j\ = \k\ shows that each central function is in terms of exactly 3r — 3 
indices, the Krull dimension of C[9ir]*^- Thus each central function 
corresponds to an admissible weight A G N^''"'^. 

Let the Clebsch-Gordan injection be denoted by 

Also, let {c*} be a basis for V^*^| ^^j^-^ and {d^} is a basis for V^(^|i|_2|A:|) 

(assuming |j| = \k\). 
Define 

■^P= (4(c:)((Xi,...,x,,)-.f(d,)))^ 

A4l'^ is a (^\i\ — 2\k\ + 1 j x (^\i\ — 2\k\ + matrix with noted s,t en- 
tries. 

In these terms we can see that 

xf(Xi,...,X,.)=tr(^f). 

Since these injections are given by iteratively using the injections 
from the rank 2 case (that is decomposing a product of tensors two 
at a time), our computation of these injections in |LP09] determine all 
such injections in general (up to a choice of associativity). 
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With this in mind, these functions take natural diagrammatic form. 
Beginning with (13.21) and its vertical reflection (providing the decompo- 
sition of the dual), tensorial contraction corresponds to gluing copies of 
the matrix variables X; in between the two diagrams. Taking the trace 
corresponds to adding a closing loop to the diagram. The resulting 
diagram is 



(3.3) xf(Xi,...,X.)=x^^,^ = 




where mi = h -\ h k+i - 2(ji H \- ji) and pi = h ^ h ii+i - 

2{ki + ■ ■ ■ + ki). The requirement |j| = \k\ becomes rrir-i = Pr~i- 

3.3. Example r = 1. The diagram is a single loop: 

The trivial representation Vq gives Xo = 1. The standard represen- 
tation Vi has diagonal matrix coefficients Xu and 0:22, hence 

Xi = Xii + X22 = tr(X). 

The remaining functions may be computed directly, or by using the 
following product formula: 

(3.4) XaXb= 

ce \a,b\ 

Explicitly, the particular case 6 = 1 is (for a > 1) 

(3.5) XaXi = Xa+l + Xa-l, 

from which the recurrence Xa+i = tr(X)Xa— Xa+i can be derived. These 
polynomials, shown in Table [H are closely related to the Chebyshev 
polynomials of the second kind. Note that the ring structure is not 
that of C[tr(X)]. 

3.4. Example r = 2. The diagram is: 

= 

Recall the decomposition 

C[SL(2,C) X SL(2,C)]^^(^'^) ^ J2 Cx;,, 

a.feeN 
eg \a,b} 
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Xo 


= 1 


Xi 




X2 




X3 




X4 




X5 


^x' 



Table 1. Rank 



1 

2x 

3x^ + 1 
Ax^ + 3x. 

Central Functions. 



where X°^ corresponds to the image of 

c c 
fe=0 k=0 

under the injection V;* ®Vc^ V* O I4* O I4 O H- 

This inclusion is determined by the Clebsch-Gordan injection l : 
Vc ^ K <8) Vft. Hence, an explicit formula for l provides a means to 
compute XI ^ directly. 

Since the general injections are determined by the rank 2 injections, 
we now review their construction. 

A few simple examples will motivate the construction of l. 

For A; = 1, 2, let Xfe = (x^j) be 2 x 2 generic matrices, and let 

X = tr(Xi) = x\i + X22, 
y = tr(X2) = xj^ +XI2, 

Z — tr(XiX2 ) = (2^11^22 "I" •^22'^ll) ~ (•^12'^21 "I" •^21'^12)- 

The map U : Vq ^ Vi ®Vi given by 

Co I— > ao (8) bi — ai (8) bo 

is invariant. 

More generally, the injection Vq ® Va is given by 

a 

U«:co^^(-ir(;)a„_„®b„. 

m=0 
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Hence, XJ] „ = ^ and ^ may be computed by: 



?> Co 

® bt - 
® ao) 
-(a? 



I* ® bo) ® (ao ® bi - ai ® bo) 
(b*®b)-(a*®a)®(b*®bi) 
ai)®(bo®b) + (at®a)®(bo 



^11 



•^22 -^12 ® •'^21 -^21 ® -^12 "I" •^22 ® -^ll 



1-^ (^11^22 •^22-^ll) (-^12*^21 "I" •^21'^12) — ^■ 

The representation ni^-y be identified with a subset of V^'^ via the 
equivariant maps 

Sym 



Proj 



where Proj o Sym = id. 

Thus, when c — a + b, l is given by the commutative diagram 



Sym 
Vr 



Va 



Proj(g)Proj 



In particular. 



0<i<a 
0<j<b 
i+j=k 



For example, consider X\ q. In this case, Co i— > ao®bo and Ci ^ ai®bo. 
Hence, 

Kfi 1-^ Co (g) Co + Ci (g) ci (ao (8) ao) (bo bo) + (a^ (g) ai) (g) (bo bo) 

1—^ x\i (g 1 + <g 1 I— ^ a^ii + ^22 = X. 



A similar computation shows that X\^^ y. 

Let 7 = l/2(a + 6-c), a = 1/2(6 + c- a), and {3= l/2{a-h + c). 
The general form of i is determined by combining these cases in the 
following diagram: 



idigiUTfglid 



Va ® Vb ^ Vp^V^^V^® Va 
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It follows that the mapping t : 14 — > K ® H is explicitly given by: 

(Dc.-- E (•)^^®[(-irc:)^7-™®b„]®(;)b, 

0<i</3 
0<j<a 
0<m<7 
i+j=k 

0<i<P 
0<j<a 
0<m<7 
i+j=k 

In |LP09] it is shown 
Theorem 3.4. Provided a > 1 and c > 1, we can write 

^a,6 ■ ^6,c-l 4a(a-l) ^^0-2,6 4c(c-l) '^o.b 

_ (a+b+c)^(a--fe+c-2)^ y.c-2 

16a(a-l)c(c-l) "^0-2,6- 

r/ie relation still holds for a = 1 or c = 1, provided we exclude the 
terms with a — 1 or c — 1 in the denominator. 

Also, note that formulae for multiplication by y and z may be ob- 
tained by applying the following symmetry relation. 

Suppose a central function is expressed as a polynomial P in the vari- 
ables X = tr(Xi), y = tr(X2), and z = tr(XiX2 so that Pa^tAv^ ^5 ^) = 
X^^(Xi,X2) for some admissible triple (a,6, c). 

Theorem 3.5. For any permutation a, 

P.(a,t,c)iy, X, z) = P„,6,,(cr"^(?/, X, z)). 

Using this symmetry and the above recursion, the ring structure is 
completely determined. As stated in the introduction, for r > 2, this 
ring structure is not known and for r = 1 , 2 it was worked out in |LP09j . 

Section [5] explores the r = 3 case using computations made with 
Mathematica. We use both the tensorial contraction method discussed 
above (which reflects our definition of central functions), and a purely 
combinatorial method that uses spin network techniques. The next 
few sections lay the groundwork for the combinatorial method, which 
comes from the representation of the central functions as spin networks. 

4. Trace Diagram Recurrences 



4.1. Gluing Lemmas. Spin networks satisfy certain "recoupling" iden- 
tities. 
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Definition 4.1 (6j-symbols, Definition 5.4 in |LP09j ). The 6j-symbols 
are the coefficients in the following change- of -basis equation: 

y - E : ^ { y ■ 

Note that this definition differs shghtly from |CFS95l IKauQl] . 
Together with the fusion identity, this provides an identity for gluing 
a strand "across" a vertex. 



Proposition 4.2 (Strand- Vertex Gluing). 

llrt ee\a.d\ ,/\ eGla.dl 



e£la,d\ if\n ee\a,d\ 

" fe\c,d\n\b,e\ 



d c f^ 




Proof. The first step applies the fusion identity (Proposition 12. 7p to the 
lower part of the diagram, introducing a summation over e. The second 
step uses the change-of-basis formula (Definition 14. ip on the upper half 
of the diagram. □ 

We now introduce specialized notation for the case of the above 
lemma with c? = 1, since it will simplify the expression of the product 
in Theorem 14. 7[ 

Definition 4.3. Given an admissible triple {a, 6, c}, a' G and 
d G [l,cj, define the fusion coefficient and the normalized fusion 
coefficient ^ by 

It is immediate that 

(4. 1) ( i i ) = v/fa'(i,a)feKi,c)^6 ( i i ) . 

Also, a' G [1, aj is equivalent to requiring a' = a ± 1, so given {a, 6, c} 
there are four choices for {a',c'}. Table [2] shows the values of the 6j- 
symbol, the fusion coefficient, and the reduced coefficient in each of the 
four cases. One can show that for any (a', c') G [1, a\ x [1, cj that 

^^{ii) = WA) and ^^,{ii)=WA)- 
These values are taken from Corollary 5.5 in [LP09j . together with 
the fact that the additional fusion constant is either fa+i(l,a) = 1 
or fa_i(l,a) = This table corrects sign errors in the 6j-symbol 

formulas provided on pp. 79-80 of [PetOGj . 
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{a', 


c'} 


lb a' a i 




^b 


{a 


+ 1, 


c+ 1} 


1 


1 


1 


{a 


-1, 


c+1} 


ec{a,b) 


ec(a,f>) 


tc(a,b) 


a 


a+1 




{a 


+ 1, 


c-1} 


ea(c,6) 


ea(c,fe) 


Ca(c,f)) 


c+1 


c+1 


^c(c+l) 


{a 


-1, 


0-1} 


e(,(a,c)(e(a,fe,c)+l) 


c)(e(a,6,c)+l) 


ei,(a,c)(c(a,f),c)+l) 


a(c+l) 


(a+l)(c+l) 


i/a(a+l)c(c+l) 



Table 2. The four basic fusion coefficients, with corre- 
sponding 6j-symbols and normahzed values. 



With a different orientation, Proposition 14.21 becomes 
Proposition 4.4. 

(4.2) Y = E ^^(a^) '^V^' • 

c'Gri,cJnrfe,a'J 

Proof. In terms of the fusion coefficient. Proposition 14.21 is 
l,.c jh Iv yc ,6 

(4.3) Y = E v(a,i)5.("^;; 

c'Grc,ijnr6,a'j 

Reflect this relation vertically using Proposition 12.31 and extend the 
strands labeled by 1 and c on both sides of the equation to obtain: 

(4-4) '((X, = E ^Aa.mbiii) " 

a'gfa.lj 
c'Gfc.lJnrfe.a'J 





By Proposition 12.51 the relation can be straightened, with the intro- 
duction of signs Si and 5*2, to: 

(4.5) ^1 Y = ^ S,BAa.m{ii) . 

^ a'era,lj ^ 

c'Grc,lJnr6,a'J 

To calculate 5*1 5*2, we must count the number of kinks in each dia- 
gram. No strands in (14.51) are kinked. The strands between h and c on 
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the left-hand side of (14.41) are kinked, giving a factor 5*1 = Sa{b, c). On 
the right-hand side, strands between b and either 1 or c are kinked, and 
so are those between 1 and c, producing a factor ^2 = Sa'{b, c')Sc'{l, c). 
Hence 



S1S2 



'l,c)Saib,c')Sa{b,c) 



V(l,a)- 



□ 



Therefore, the coefficient is (v(l, a)) db ( i i)=dbi i i)- 

The reflective symmetries guaranteed by Proposition 12.31 imply that 
corresponding formulas for the products 

V/' A 

can be generated by reflecting (14.31) . In particular, the terms in these 
summations have the same coefficients. 



4.2. Simple Loop Recurrences. Relation (14.31) can be "stacked" to 
obtain a more general formula for gluing across two or more vertices. 

Lemma 1. // the edges labeled by {ao, . . . , a„} are disjoint, then 



ry {a',eri,adnr6.x-iJ}Li 



f,j,(ao,l)f<(a„,l)n5au(«-i>l)^''. ( ab i 



i=l 




Proof. For clarity, we present the concrete n = 2 case here, from which 
the general pattern can be seen. First, apply (14. 3 p on the upper and 
lower halves of the diagram separately: 




aoefi.aoj 

a'^gfl.aijnrbi.a^j 
a2€\l,a2\n\b2,a\\ 



E 



{<Gri,adnrb»,<_iJ}? 




1' 
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Note that the bubble identity (12.101) allows us to use the same index 
a\ above and below the bubble in the second diagram and contributes 
the ba'^(l,ai) term in the second step. Using (14. ip and the fact that 
fc(a, 6)bc(ci, &) = 1, the coefficient can be expressed as 

2 

^f,,(ao,l)f4(a2,l)n5au(«-i'l)^^. i) ■ 

i=l 

The case n > 2 follows by induction. □ 

We now introduce notation that describes relabelling of spin net- 
works. In what follows, we let s be a closed trivalent graph with edges 
E = (ei, 62, ... , Cm)- Given a labelling I : E N of the edges by natu- 
ral numbers such that the labels at all vertices form admissible triples, 
let Si denote the resulting spin network. 

Definition 4.5. An admissible relabelling of I is a labelling la of a 
subset of edges a : (cj^, Cj^, ■ ■ ■ , ej„) N such that the resulting network 
Si^ is still admissible. 

The next theorem says that multiplication of a spin network by a 
simple loop may be expanded in terms of diagrams with admissible 
relabellings. 

Definition 4.6. A simple cycle (ei, 62, . . . , e„) is an ordered tuple of 
distinct edges in a graph that begins and ends at the same vertex. 

Theorem 4.7 (Simple Loop Multiplication Formula). Let s be a triva- 
lent graph with admissible labelling I : (ci, 62, . . . , e^) — > N. Let 7 
denote both a simple cycle (e^j, ej^, . . . , Cj^) in s and the corresponding 
spin network. Define ai = l{ej.), and let hi G [aj_i,ajj be the third 
label on the vertex joining edges ej._j and Cj.. Let a' : cj^ ^-^ a- be an 
admissible relabelling of si . With the understanding that Cj^ = cj^ and 

{<eri.«>Jnrb>,a;_iJ}?=i Vi=i / 

where Ni is the sum of (i) the number of times 7 crosses the edge ej^, 
(a) the number of local extrema along the edge Cj. that do not occur at 
a vertex, and (Hi) the number of times Cj. adjoins a vertex as shown in 
Figure \^ 
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1, .a 



8 

(a) (b) 



7iv«i ibi±i bi±i. ai -y 7, ,ai±i 



(c) 



Figure 4. Local contributors to sign of simple loop mul- 
tiplication formula. 

Proof. First, suppose the loop has the form of Lemma [H with the oq 
and a„ edges coinciding, as follows: 



Popping the final bubble introduced along these edges results in an ad- 



ditional factor of ba(,(oo, 1), which cancels with y'fa(,(ao; l)fa^('2n, 1) = 

fajj(ao, 1); the last equality holding since a„ = a^. So in this case the 
coefficients of the summation are 



-1 



i=l 

In general, the relative positions of the labels {ai-i,ai,bi} in the 
vicinity of a vertex matters. If the product appears locally at a local 
extrema, as in (14. 2 p or its reflection, then the term Sa'.{l,ai) is un- 
necessary. Otherwise, in the cases depicted in Figure |l](c), the term 
remains. 

If a crossing or local extrema occurs along an edge, then the "bubble 
popping" step in Lemma [T] becomes one of the following: 




Sa'(l,ai)ba'(l,ai] 



ai 



2 1 

1 



aiy 1, .1 ,ai 



Sa'(l,ai)ba'(l,ai) 




The signs are calculated using Lemma [2741 and the stronger Proposition 
2.5[ In the second case, the sign is calculated by comparing the ta'- (1, Oj) 

kinks in the diagram VSJ I with the diagram 'I y , which has 

no kinks. In each case the additional sign Sa^(l,aj) adds one to the 
exponent iVj. □ 
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Rearranging the terms in the above theorem provides a recurrence 
formula, in which each diagram can be written in terms of diagrams of 
lower order, where order is defined as follows. 

Definition 4.8. The order of a spin network si is the sum YllLi K^i)- 

Corollary 4.9. Let s, I, 7, Oj, bi, N.i, and a' he defined as in Theorem 
\4. 7[ If is an admissible partial relabelling, where a{ejj = l{ej.) — 1, 
then Si can be expressed in terms of'y and spin networks of lower order: 
(4.7) 

= rsu- E (fl^a'Sha. - 1)^^^, s,,. 

{a'^e\l,a,-l\n[b„ar_^]}f^^ \i=l / 
some a'^j^ai 

Proof. The term with a'(ej-) = in Theorem 14.71 is the unique one 
with highest order. So one may re-index by replacing the with — 1 
as shown. Since 3^,(1,0, -!) = +! and ^b, {aZi-i = 1' there is 
no coefficient on the term. □ 

4.2.1. Some Examples. We illustrate the application of these theorems 
in a few basic examples. Note that the sign is only a factor when 
a[ = Qi — 1 since 

X J+1 a • = aj + 1 
1^-1 a. = aj-l. 

We will give the recurrences in their most general form; for some choices 
of labels, the non-admissible terms should be excluded. 



Example 1. A single edge loop has a two-term recurrence: 




The loop 7 consists of a single edge, with Ni = 1 since there are no 
crossings and 7 has one extremum that does not occur at a vertex. 



Setting 6 = 0, one obtains a formula equivalent to fl3.5l) . which can 
be used to compute the rank one central functions. 
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Example 2. A two-edge loop has the following four-term recurrence: 

61^ 



(4.9) 



0.1/ 



a[=ai±l \i=l 
a2=a2±l 



61' 



Since A^i = + + = A^2; no additional signs are necessary. 

4.2.2. Application to Rank Two Central Functions. The rank two cen- 
tral function is 

^:,.(Xi,x2) = 

There are three simple loops: (a, b) corresponding to tr(XiX2 ^), (a, c) 
corresponding to tr(Xi), and (6, c) corresponding to tr(X2). Each loop 
provides a different recurrence, a fact which was used in |LP09j to 
obtain a new proof of a classical theorem of Fricke, Klein, and Vogt. 
For example, closing off each term in ([2]) gives the recurrence 




(4.10) 




n + 1, 



,6 + r 




+ 



ea(6,c) 



2 + 1, 



-1) 




+ 



et.(a,c) 
a(a+l) 



2 a — 1, 



i + r 




+ 



tcia,b)^(tia,b,c)+l)^ 
a{a+l)b{b+l) 





Note the similarity to the formula in Theorem 13.41 

4.2.3. Application to Rank Three Central Functions. The left-associative 
rank three central functions are 



Xa,b,c,d,e,f 



There are six simple loops in the diagram: 

(a,6)^tr(XiX2^), {c,d)^tT{X,), 
(a, e, d, f) ^ tr(Xi), (6, e, rf, /) ^ tr(X2), 
(a, e, c, /) ^ tr(XiX3 1), (6, e, c, /) ^ tr(X2X3-i). 

However, as indicated in the appendix, up to seven variables may be 
required in the expansion of rank three central functions, so the simple 
terms do not suffice to compute all central functions. This case will be 
treated in detail in section 15.21 
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4.2.4. Application to General Central Functions. In section [3], the left 
associative central functions were shown to be given by 




Theorem 4.10. There are lil^ii possible simple loops in the general 
rank r left- associative central function. Multiplication by any of these 
possible loops gives rise to a recurrence as described in Theorem \4.'T\ 

Proof. Simple loops must pass through precisely two of the edges la- 
beled by i;, implying that there are {^~^) = ^^^^^^i^ possible simple 
loops. □ 

The simple recurrences given in this way correspond to the following 
subset of generators of C[^r]^: {tr(Xi), ...,tr(X^)} U {tr(XiXj) | 1 < 
i < j < r} (see appendix for a full description of minimal generators). 



4.3. Non-Simple Loop Recurrences. As mentioned in the previous 
section, multiplication by simple loops does not provide sufficient re- 
currence relations to reduce any trace diagram to its simplest pieces. 
In the rank three case, for example, multiplication by tr(XiX2 ^Xs) 
requires a non-simple loop: 




While Corollary 14.91 can be applied to any simple loop in a trace dia- 
gram, it is not sufficient to compute the value of an arbitrary diagram. 
The simplest case where it fails is the "barbell" depicted in Figure [51 
There is no simple loop recurrence, since subtracting one from either 
loop produces a non-admissible diagram. However, the non-simple loop 
depicted on the right of Figure [5] does provide a recurrence. 
Define 
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Figure 5. Trace diagram with no simple recurrence 
(left). The non-simple loop which can be used to pro- 
duce a recurrence (right). 



When applied to the simple loops in the diagram, (14.81) implies 



tr(X) ■ Xl. = Xa+l,c - ( ""a' "a' ) " Xa- 



These equations can be rearranged to provide recurrences for the X 
and Y loops. The third recurrence corresponds to the multiplication 
by tr(XY): 



Proposition 4.11. 




-Y: v(i, a)v(i, c)L ( i 'I ) #c ( ^ i ) h ( ) dc ( i ) 



6" 



b'G\b,l\n\a,a'\ 
c'G[c,l\r\lc,b'\ 
b"£\b',l\nla,a'\nlc,c'\ 



Proof. The computation is a straightforward application of Lemma 
and the bubble identity (12.101) . For simphcity, the matrices are not 
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shown on the diagram; nothing changes if they are included. 



c'e\c,i\n\c,b'] 





= E ^4 ) ( c i ) ^')^a' ( ^; ) 5,. (i, b')^, ( ^ ^; ) (i, 

6'erf',iJnra,a'J 

c'€\c,l\nlc,b'\ 
b"e\b',l\nla,a'\n\c,c'\ 

b'€lb,l\nla,a'\ 

c'e[c,i\nlc,b'\ 

b"e\b',l\n\a,a'in\c,c'i 

□ 

The exphcit computation of these functions is shown in Table |3] with 
X = tr(X), y = tr(Y), and z = tr(XY). We omit the cases where 
6 = 0, since they can be obtained as a product of two rank one central 
functions. 



a+b+c 


-b 

A.a,c 


formula 


4 


xli 


^ 2 


5 


xli 


xz-^^ 




xh 


9 

yz-^ 


6 


xli 


™2 x'^^y _|_ xy 2z 
X ^ 2 3 3 




xh 


12 2 

xyz — i^x y 






„,2 xy^ 1 xy 2z 
y ^ 2 3 3 


7 


4 cases 




8 


X2,2 

5 more cases 


Z xyz +6+3 + 3 3 



Table 3. Barbell Functions. 
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4.3.1. The General Case. Other non-simple recurrence formulas may 
be similarly derived. In every case, the coefficients will be signed sum- 
mations over the normalized fusion coefficients, with the sign depending 
on the particular configuration of the graph. 

5. Rank Three Central Functions 

5.1. Tensorial Algorithm. We now heuristically describe an effective 
algorithm to compute central functions for ranks 1,2, and 3 using a 
"tensorial contraction" method. This algorithm has been successfully 
implemented in Mathematica, and a semi- documented "notebook" is 
available to the reader (at the first author's website). 

Let t ^ C^'- be the mapping given by specifying A^^ minimal 
generators for C[9^r]*^5 let X,. = tiJRr) be the image of t. The projection 
yir is dual to the inclusion C[9^r]'^ C C[9lr-]- Denote by [p] the 

equivalence class of p G IHr where p p' ii and only if Gp fl Gp' ^ 0, 
that is their conjugation orbit closures intersect non-trivially. can 
be understood as the space of such equivalence classes. In the case 
r = 3 this mapping is given by [p] ^ 

(tr(p(xi)), tr(p(x2)), tr(p(x3)), tr(p(xiX2)), tr(p(xiX3)), tr(p(x2X3)), tr(p(xiX2X3))) . 

In |Gol08j an explicit global slice to this projection is constructed. Us- 
ing this slice and an iterative application of the work in |LP09] we 
compute rank 3 central functions tensorially and describe this algo- 
rithm below. 

To ease the notation let = tr(Xfc) = tr(p(xfc)), tjj = tr(XjXfc) = 
tr(p(xjXj)), and tyfc = tr(XjXjXfc) = tr(p(x^jXfc)). 

The projection vr : X,. ^ given by [p] ^ (^1,^2,^3,^12,^13,^23) 
is a branched double cover. In particular, there exist polynomials in 
P, Q G C[ti, t2, ^3, ^12, ^13, ^23] SO that 

C[X3] ~ C[ti, t2, ts, tl2, tl3, ^23] [il23]/(43 " ^^123 + Q)- 

The two roots of the irreducible generator of the ideal are related by 
the formula: 

^132 = —t\t2tz + ^12^3 + ^2^13 + ^1^23 ~ ^123- 

This provides the formula for P; the formula for Q is 

^123^132 = (i^l+i^2+^3) + ('^12+^23 + '^13)~(^1^2i^l2+^2i^3'^23+^3i^l^l3)+^12i23^13~4. 

To construct a slice we must construct a triple of matrices (Xi, X2, X3) 
for every 7-tuple in the image of t. 
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Let Xi = ^ 0^1 X2 = ^ ^ 1 ^ 1 5 where w+ 1/w = ti2- 
Then letting 



where s is a solution to the equation det(X3(s)) — 1 = 0, will give the 
desired slice (see |Gol08] ). 

5.1.1. Step 1. Define tensor products and duality scalar products. The 
scalar product satisfies 

n:_fc(n„-0 = ^, Ski = ^ki/Q, 

for Hn-k e Vn. 

5.1.2. Step 2. Define an algorithm which determines the form of a lin- 
ear representation of an element in SL(2, C) on the symmetric power 

Vn. 

This comes from 

n*_;,^(X ■ rtn-i) = n*_fc ((xnCi + X2ie2)''"'(xi2ei + X22e2)') 



E a) "MO 



•^11 -^12 •^21-^22- 



0<i<n-l 
0<j<l 



For instance, if 5^ = ^ ^ ^ ^ , then the induced action on V2 is given 

ah 6^ \ 

by I ac \{hc + ad) hd , and on V3 by 
C(i d'^ I 

( a^ a% a}? \? \ 

a^c I {da'^ + 2bca) | (cfe^ + 2adb) }?d 

ac^ \{hc^ + 2adc) l{ad'^ + 2hcd) bd^ 

y c^d cd"^ d^ J 

Upon contracting tensors, we will need to know the polynomial ma- 
trix coefficients of the action on a given Vn- This routine will allow us 
to read off such matrix coefficients, since the pairing for n*_^(X ■ n^.;) 
just becomes the n — k + l,n — I + 1 entry of the symmetrized matrix. 
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5.1.3. Step 3. Given a triple (a, 6, c) determine all admissible 6-tuples, 
and mark them by an enumeration of the multiplicities. We iteratively 
decompose tensor triples K ® H ® K using a left associative iteration 
algorithm. In other words, we decompose (V^ ® VJ,) ® Vc, using the 
decomposition of ® VJ, ~ Va+b-2k- Then for each allowed value of 
< k < a+b, we decompose Va+b-2k'^Vc in the same fashion. We know 
that (a, b, a+b— 2k) is admissible if and only if \a—b\ < a+b— 2k < a+b. 

For instance Vd injects into V3^V2<S>V2 if and only if o? e {7, 5, 3, 5, 3, 1, 3, 1}. 
The multiplicity arises since 

(^3 ® V2) (8)^2- (V^ © © Vi) 0V2^V50V2®V30V2®Vi0V2 

« 14 © 1/5 © © ■^4; © ■^^3 © © © Vi. 

5.1.4. Step 4- Define the injection of basic elements. This is done using 
the formula for the mapping t : Vc — >■ © H- 

0<j</3 
0<i<a 
0<m<7 
i+j=k 

0<i</3 
0<j<a 
0<m<7 
i+j=k 

and iterating for each summand. 

Also since we are using left associations, i.e. grouping Va and V5 
in Vi © Vf, © we first inject Vd ^ Va+b-2k © K and then inject 
Va+b-2k © V; ^ K © 14 © V; using the injection Va+b-2k "-^ K © H 
tensored with the identity. 

Likewise, we have an injection for the dual injection ^ V* © 

5.1.5. Step 5. The central function begins in End(Vd)^'"*^^''^) and then 
is mapped to {V^ © Vaf^^'^''^^ by 

d c 
fe=0 fe=0 

The central function is then determined by the composite injections 

End(Vd)S^(2'^) ^ {V^^Vdf''^'''^^ {V:®V,*®V:®Va®Vb®V,f^^''^'^\ 

Once we realize the exphcit form of the injections ^ K © H © K 
we can write the central function $^^=0 ^ ^ central tensor in 
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terms of Va Vb ® Vc- We note that the coefficients are a bit dehcate 
here: ffist one chooses to include on only one the vectors dk or 
d^, but not both, for each summand of the central function. Second 
one must make sure to include a 1/ ) when including each factor of 
each summand (both the vector and its dual) into a tensor Va (S> H. 
We observe that the rank 2 coefficients cancel, but for rank 3 they 
generally do not cancel since we are iterating injections. Lastly, we 
have included the dual pairing binomial coefficients in our expression 
of the symmetrization of a generic matrix, so we do not need to include 
any further binomial coefficients (there are a lot included already!). 

5.1.6. Step 6. Once we have the central function as a linear combina- 
tion of tensors, we must "contract the tensor" using the mapping that 
associates a linear combination of tensor (of the type we have been con- 
sidering) to a polynomial function in the coordinate ring C[!Hr.]. Recall 
that this is determined by mapping (e*^ ® e*^® ■ ■ ■ ® e*^) ® {ej^ ® ej^ ® 
■ ■ ■ <S) Cj^) to the polynomial function 

(Xi,X2, ...,X^) ^^ e*^(Xi ■ ejJe*2(X2 ■ e^J ■ ■ -e^X^ ■ e^J- 

Once a generic matrix Xj is mapped to an automorphism of Vd, call 
it Sym'^(Xj). Then e*^(Xi ■ e^J is just the (ii + l,ji + 1) entry of 
Sym'^(X,). 

5.1.7. Step 7. Lastly, using the Goldman slice (see |Gol08] ) 

Xs = SL(2,C)^V/SL(2,C) ^^3 = SL(2,C)^^ 

to the categorical projection 9^3 — > X3, we express the invariant poly- 
nomials (which are polynomials in the generic matrix entries) from the 
tensorial contraction in terms of the seven invariants: 

{tr(Xi), tr(X2), tr(X3), tr(XiX2), tr(XiX3), tr(X2X3), tr(XiX2X3)}. 

It is only this last step that does not generalize to r > 4, since we 
do not have a slice. However, for r > 4, we can still contract the 
tensor (Steps 1- 6) and then use a Grobner Basis algorithm to express 
the polynomials in terms of traces (minimally in fact since minimal 
generators are known for all r > 1). 

5.2. Combinatorial Algorithm. Trace diagrams allow a more com- 
binatorial approach to computing the central functions. The key point 
is that any central function can be reduced in terms of simpler central 
functions. As discussed in section 14. 2[ a recurrence formula exists for 
each loop in the diagrammatic depiction of the central function. 
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In this section, we show how this process works for r = 3 diagrams 
of the form: 



X,a,b,c,d,e,f 

We call the sum of a + 6 + c the order. Then the algorithm reduces 
the diagram into a sum over diagrams in a lower order. The algorithm 
terminates with the base case Xo,o,o,o,o,o = 1- The process uses eight 
different recurrence formulas, corresponding to the following trace vari- 
ables: 

{tr(Xi), tr(X2), tr(X3), tr(XiX2-i), tr(XiX3 1), tr(X2X3-^), 

tr(XiX~^X3),tr(XiX3X-^)}. 

The final two recurrences, which correspond to nonsimple loops in the 
diagram, require particular attention. 

5.2.1. Algorithm Overview. As mentioned in section the recurrence 
corresponding to a loop with n edges contains a maximum of 2" terms. 
For that reason, it is more efficient to begin the computation with the 
shortest loops in the figure. The algorithm implemented in Mathemat- 
ica begins with the following steps: 

(1) If te{a,b) > and Cf{a,b) > 0, reduce along the {a,b) loop, 
corresponding to tr(XiX^^); 

(2) otherwise, if Ce(c, d) > and C/(c, d) > 0, reduce along the (c, d) 
loop, corresponding to tr(X3); 

(3) otherwise, if Cf,(a, e) > 0, efe(a, /) > 0, tc{d, e) > 0, and tc{d, f) > 
0, reduce along the (a, e, d, f) loop, corresponding to tr(Xi); 

(4) otherwise, if ta{b, e) > 0, ea(6, /) > 0, tdd, e) > 0, and tdd, f) > 
0, reduce along the (6, e, d, f) loop, corresponding to tr(X2); 

(5) otherwise, if Ca{b, e) > 0, Ca{b, f) > 0, td{c, e) > 0, and td{c, f) > 
0, reduce along the {b, e, c, /) loop, corresponding to tr(X2X3"^); 

(6) otherwise, if eb(a, e) > 0, Cb{a, f) > 0, td{c, e) > 0, and erf(c, /) > 
0, reduce along the (a, e, c, /) loop, corresponding to tr(XiX3"^); 

(7) otherwise, if e = 0, reduce along the nonsimple (a, b, e, c, d, e) 
loop, corresponding to tr(XiX3X2^); 

(8) otherwise, if / = 0, reduce along the nonsimple (a, b, /, c, d, f) 
loop, corresponding to tr(XiX2^X3). 

Proposition 5.1. Any admissible Xa,b,c,d,e,f can be reduced via one of 
the above cases. 
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Proof. Consider the following set of edge types: 

(5.1) {Ca{b, e), tb{a, e), ta{b, /), tb{a, /)}. 

If any two of these is zero, then the central function is reducible. To see 
this, suppose without loss of generality that e-a{b, e) = 0. If efe(a, e) = 
then e = ta{b,e) + tb{a,e) = 0. If tb{a, f) = 0, then a = b + e and 
6 = a + / so that e = and / = 0. Finally, if Ca{b, /) = then 
a = b + e = b + f and so te{a,b) = 2b = tf{a,b). So either 6 > 
and reduction by the (a, b) loop is possible, or 6 = 0. In this case, the 
diagram is reducible by either tr(Xi), tr(XiXj^), or tr(X3). 
The same logic can be applied to the set of edges 

(5.2) {Cc{d, e), td{c, e), tc{d, /), td{c, /)}. 

But note that if at most one of flS.ip and at most one of fl5.2l) is zero, 
then one of the four reductions corresponding to a loop of length four 
is possible. □ 

5.2.2. Reduction along Simple Loops. The simple loop reductions may 
have up to 4 terms in the tr(X3) and tr(XiX^^) cases, and up to 16 
terms in the other four cases. While this is a large number of terms, 
writing out the recurrence is a straightforward process. Each is an 
immediate consequence of Theorem 14.71 For example: 



tr(XiX2 '^)Xa,b,c,d,e,f — Xa+l,b+l,c,d,eJ 

, f u\ ( u\^a{b,e)ta{bj) 

+ Se(a, b)5f[a, b) _^ Xa+l,fe-l,c,d,e,/ 

, ( n ( u\^b{a,e)tb{a, f) 

+ Se[a, b)Sf{a, b) Xa-l,fe+l,c,d,e,/ 

ee(Q, &)e/(a, &)(e(a, &, e) + l)(e(a, &, /) + 1) 
a[a + l)b[b + l) 

The remainder of the recurrences may be found in a Mathematica note- 
book written by the second author (available on his website). 



5.2.3. Reduction along Non-Simple Loops. In the final cases, either e = 
or / = 0, and the central function is topologically equivalent to a 
barbell. Proposition 14.111 could be adapted to this case, but we prefer 
to show directly how the rank three central functions are related to the 
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barbell functions Xa c- When / = 0, the diagram is 




The sign (— l)^ arises from Proposition 12.51 since there are | kinks in 
the left diagram and none on the right. So x{(^: b, c, d, e, 0) may be 
obtained from Table [3] using the substitutions 

X ^ tr(X3) = ts, 

y ^ tr(X2 'Xi) = ht2 - ti2, 

Z (-^ tr(X3X2 ^Xi) = ti23 + tit2t3 — ti2t^ — tlt23- 

Similarly, when e = 0, the diagram is 




The proper substitutions in this case are 

X ^ tr(XiX2 ^) = tit2 - ti2, 
y t-> tr(X3) = tg, 

Z ^ tr(X3XiX2 = h3t2 - tl23. 



5.3. Computations. Recall that every rank 3 central function has 6 
indices. The first 4 come from the inclusion of a symmetric tensor 
Vd ^ Va ® Vb <^ Vc- Suppose it occurs with multiplicity m. Then for 
^ ^ h j ^ fn, the central functions are indexed by X^abcd^ ^^^^ ^ 
j for V^. Any 6-tuple arising in this fashion is called admissible. Note 
that this choice of index differs from Xa,b,c,d,e,f- By letting c = we 
recover the rank 2 central functions in |LP09j and by letting b = c = 
we recover the classical rank 1 central functions. 

We call the set of all admissible 6-tuples satisfying a + b + c = s 
the s- order. Table H] shows the first four orders of rank three central 
functions. 
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0- order (1 function) 

1,1 _ 1 
Xo,o,o,o — ^ 

1- order (3 functions) 

1,1 _ , 1,1 _ , 1,1 _ , 
Xi,o,o,i — f^i Xo,i,o,i — Xo,o,i,i — ^^3 

2- order (9 functions) 

V-l'l _ +2 _ 1 ^1,1 _ hts I 1,1 _ t2t3 I t23 

A2,0,0,2 ~ ''1 ^ -^1,0,1,2 ~ 2 2 ^^0,1, 1,2 ~ 2 2 

Xi;i,0,2 = ^ + ^ XlAl.O = ^1^3 - ^13 XoXl,0 = ^2^3 " ^23 

Xl4,0,0 = ^1^2 - ^12 Xo'2,0,2 = *2 ~ 1 Xo!o,2,2 = ^3 ~ ^ 

3- order (20 functions) 

X3,0,0,3 = ^1 ~ 2ti Xm,1,1 = — 1^1^2^3 + + ^1^23 — ^123 

,,1,1 Ij- J-2 I 2ti2f 1 2*2 A/1,1 3 J. J. J. I tl2t3 fatia ^1*23 

X2,l,0,3 - 3''2'^1 ^ 3 T Xl,l,l,l - 4111213 + 2 2~ 

^1,1 — + +2 ^ ^ t2 ^1,1 _ U +2 I 2ti3t3 2ti 

X2,i,o,i — '^2ti - 11211 - Y Xi,o,2,3 — a'-i^s 3 3" 

A.1,1 _ ly- v-2 I 2ti3ti 2*3 ,,1,1 — + +2 y. + tl 

X2,0,l,3 — s^'S'^l 3 3~ Xl,0,2,l " ^1^3 ~ ^^1313 " 2" 

X2io,l,l = ^3^? - ^13^1 - y Xo^3,0,3 = ^2 ^ 2t2 

^1,1 _ If f2 I 2*12*2 _ 2*1 1,1 _ U ^2 I 2*23*2 _ 2*3 

Xl,2,0,3 ~ 3''1''2 I 3 3 A0,2,l,3 ~ 3''3''2 3 3 

Xl'2,0,1 = ^1^2 ~ ^12^2 — Y Xoi2,l,l = ^3^2 ~ ^23^2 ^ y 

,,1,1 _ f3fl2 _|_ *2*13 _|_ *1*23 ,1,1 — ly- +2 I 2*23*3 _ 2*2 

Al, 1,1,3 ~ 3 "I" 3 "I" 3 A0,l,2,3 ~ 3''2''3 3 3 

Xl4,l,l = ^1^2^3 — ^3^12 Xo!l,2,l — ^2*3 — ^23*3 — ^ 

XlXl,l ^ 2^1*2^3 - - ^2^13 + ^123 XoA3,3 = ^3 ~ 

Table 4. Order of rank 3 central functions. 



A more complicated example is 

21 "^23 ^223 "^23 ^ 3 2^j^ ^22 

X3i2,2,3 = ^^^2^1 + J^^2^3^1 ~ ^^^3^1 ~ g^2t3^23^1 + ~ Y5^2^3^12^1 

^22'^ 2'^ 2^ 2^ 2 

~ T7T^2^3^13^l + 77^3^13*^1 + 7'!^3'^12^23^l + T;T^2^13^23^i + 77^2^3^123^1 

10 15 5 10 3 

1 1 ? 1 9 1 

2 2 22 2 22 2 

- -^23^123^1 - -^2^1 + Y^^2^3^1 + Y5^3^1 ~ ^^^3^12^1 ~ 30^^^^^ 

^22 "^^2 ^^2 

+ g^2^13^1 ~ 15^13^^ ~'~ 30^^'^^^ ~ g^2^3^12^13^1 ~ 2^2^3^23^1 

1 1 1 ti 4 2 
+ 7^77^12^13^23^1 + 77^3^12^123^1 ~ 77^2^13^123^1 ~ — + 77^2^12^13 

30 3 z io 10 

+ 15^'^'^'' " 1^ + 3^^^^'^^^^ + Y5^''^'^'' -^-W' m^'^''^'' 

2 11 ^23^123 
+ 77^2^13^23 ~ 77^2^3^123 ~ 7^12^13^123 H 7; ■ 

15 3 6 6 
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Using both our tensorial algorithm and our combinatorial algorithm 
we have computed all orders up to 10 which gives 2254 known exam- 
ples. We note that on the same computer it took 2 minutes to compute 
all rank 3 central functions up to order 10 with the combinatorial al- 
gorithm, but took in excess of 6 hours to do so with the tensorial 
algorithm. 
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Appendix A. Background 

A.l. Some SL(2,C) Representation Theory. We now review some 
basic SL(2, C) representation theory, and reintroduce some notation 
from our first paper [LP09j . 

Let Vq = C = Vq he the trivial representation of SL(2,C). Denote 
the standard basis for by ei = ( J ) and 62 = ( ? ) , and the dual basis 
by its transpose: el and 62- 

Then the standard representation and its dual are 

V = Vi = Cei® Ce2 and V* = V{ = Cel © Ce^. 

Denote the symmetric powers of these representations by 

Vn = Sym"(y) and V* = Sym"(V*). 

One can show Vn ~ (^)* ~ SL(2, C)-modules, which is particu- 

lar to SL(2, C) and not obvious. 

Proposition A.l. The symmetric powers of the standard representa- 
tion of SL(2, C) are all irreducible representations and moreover they 
comprise a complete list. 

For proof see [FH91j . 

The tensor product Va ® H, where a, 6 G N, is also a representation 
of SL(2, C) and decomposes into irreducible representations as follows: 

Proposition A. 2 (Clebsch-Gordan formula). 

min(a,b) 

(A.l) Va+b-2j. 

j=0 
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The particular Va+b-2j summands in this formula are described by 
the following: 

Definition A. 3. Given a,b G N, we write c G \a, b\ and say that 
{a, b, c} is an admissible triple for all c = a + b — 2j , < j < min(a, b). 

We remind ourselves of Schur's Lemma for later use: 

Proposition A. 4 (Schur's Lemma). Let G be a group, V and W irre- 
ducible representations of G, and f G HomdV, W) with f ^ 0. IfV^ 
W, then dime HomG(V^, W) = 1; and ifV ^ W , then dime HomG(V", W) 
0. 

A tensor product vi ® V2 ^ ■ ■ ■ ® Vn G V^"' projects to Vn by sym- 
metrizing. We define its image under this operation by 

where the sum is over all permutations on n elements. There exist 
bases for Vn and V*, given by the elements 

On-k = e^~''e^ = ci o d o ■ ■ ■ o ci o Co o Co o ■ ■ ■ o Co and 

n—k k 

^n-k = i(ii) (ea) =eioei0...oei 0630 620... 62, 



n—k k 



respectively, where < k < n. These elements are described in dia- 
grammatic form in section I2.1[ 

The "dual" pairing between Vn and V* is given by 

where S^i is the symmetric group on n elements. In particular, 

(n — kV.kl ^ e / 

<-ki"n-l) = ^, = 5kl/Q. 

Let g = (^^^ G SL(2,C). The SL(2, C)-action on K is given 

\5'21 ^722 / 

by 

g ■ nn-k = {guei + g2ie2T''' {guei + ^'2262)'' 



E (?) (9u'-'9'l2-W2i922) n 



0<j<n-k 
0<i<k 



1n-(i+j)- 



36 SEAN LAWTON AND ELISHA PETERSON 

For the dual, SL(2, C) acts on V* in the usual way: 

A. 2. SL (2, C)-character varieties. Let C[x*j]/ A be the complex poly- 
nomial ring in 4r variables {1 < k < r and 1 < i,j < 2), where A is 
the ideal generated by the r irreducible polynomials 

It is not hard to see that C[9^,.] = C[xf^]/A. Let x^^ be the image of 
x'^j under C[x^j] C[xfj]/A. Define 

Xi2 
■^21 -^22 

Such elements are called generic unimodular matrices. We note that 

■^11 -^12 
■^21 -^22 

are simply called generic matrices. Let G = SL(2, C). G acts on C[lHr] 
as follows: 

g ■ xlj = where := gX^g'^, and g e G. 

The ring of invariants CfD^r]*^ is a finitely-generated domain (see Nag64| ) , 
which implies 

Xr := Spec_, {C[y\rf) 
is an affine variety over C, called the G-character variety of F^. It is 
the variety whose coordinate ring is the ring of invariants, that is, 

C[Xr] = C[^.f. 

Closely related to C[Xr] is the ring of invariants 

C[2),] := C[0[(2,C)^n'^'^^''^^ = C[4f L(2,c)_ 
In fact one can show (see |Law08] ) 

C[2),]/A ^ C[X,] 

Otherwise stated, 

C[4.]SL(2,c)/^^(C[xJ]/A)^^^^'^\ 

which is true because SL(2, C) is linearly reductive. In 1976 Procesi 
proved (in the context of n x n generic matrices) 

Theorem A. 5 (Procesi). C[2).r] is generated by the invariants 

{tr(Xj^Xj2 ■ ■ -Xjj^)}, 
where Xj are generic matrices. 
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Evidently, this ring is multigraded. Finding minimal generators 
amounts to finding all linear relations among generators of the same 
multidegree in the vector space 

v. = c[2).]+/(q?)r]+)', 

where C[2)]''' is the ideal of positive terms. It is worth observing that 
there can be no relation among generators of differing multidegree. 
Any such linear relation in Vr pulls back to a reduction relation in 
C[2)r] which in turns projects to a "reduction relation" in C[Xr]. The 
Cayley-Hamilton equation gives 

-tr(X)X + det(X)I = 0. 

And if we assume det(X) = 1, as is the case in C[Xr], we easily derive 
tr(X-^) = tr(X) and tr(X2) = tr(X)2-2. Hence the generators tr(X2) 
in C[2)r] project to tr(X)^ — 2 in C[Xr] and so are freely eliminated. 

We work a couple examples before moving on: C[2)i] is algebraically 
generated by {tr(X")}. But if ri > 3 the Cayley-Hamilton equation 
provides relations which express the generator in terms of the gen- 
erators tr(X) and tr(X^). However, the dimension of this variety is 
computed to be 2. Thus there can be no further relations. 

Let Ai{R) be a set of minimal generators of a ring R. Then in 
[LawOSj it is shown, as with a generalization to SL{n, C), that 

M{Cm) - {tr(X?), tr(X^)} = M{C[Xr]), 

as long as the generators of C[2)r] are taken to be of the form in Pro- 
cesi's theorem. 

Multiplying the Cayley-Hamilton equation on both sides by words U 
and V allows us to freely eliminate the generators of type: tr(UX'"V) 
as long as n > 2 and at least one of U or V is not the identity. So for 
the case C[2)2], we are left with the generators 

{tr(Xi), tr(X2), tr(X2), tr(X2), tr(XiX2)} 

since any other expression in two letters would result in a sub-expression 
with an exponent greater than one, which we just showed was impossi- 
ble. Since in this case the dimension of the variety is 5, there can be no 
further relations and thus these generators are minimal and 2}r ~ C^. 
We can conclude that Xi ~ C and X2 ~ C'^. 

More generally, it can be shown that there are no generators neces- 
sary that have word length 4 or more (see |Gol08j for an exposition). 
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In particular, 

tr(WiW2W3W4) =tr(Wi)tr(W2)tr(W3)tr(W4) + tr(Wi)tr(W2W3W4) 
+ tr(W2)tr(W3W4Wi) + tr(W3)tr(W4WiW2) 
+ tr(W4)tr(WiW2W3) + tr(WiW2)tr(W3W4) 
+ tr(WiW4)tr(W2W3) + tr(WiW2)tr(W3W4) 

- tr(WiW3)tr(W2W4) - tr(Wi)tr(W2)tr(W3W4) 

- tr(Wi)tr(W4)tr(W2W3) - tr(W2)tr(W3)tr(WiW4) 
(A.2) -tr(W3)tr(W4)tr(WiW2). 

and 

tr(XYZ) + tr(XZY) = tr(XY)tr(Z) + tr(XZ)tr(Y) 
(A.3) + tr(ZY)tr(X) - tr(X)tr(Y)tr(Z). 

Equations flA.2l) and flA.3l) together imply there are Nr = ^ g^^^ 
minimal generators of C[Xr]- They are: Qi = {tr(Xi), tr(Xr)} of 
order r, Q2 = {tr(XjXj) | 1 < « < j < r} of order and Q3 = 

{tr(XiX,Xfc) I 1 < z < J < A; < r} of order ±zlMlzll_ 

Enumerating these generators, {ti, ...,tiy^} defines a polynomial 
mapping 

t = (ti,...,t^J:X, ^C^'-. 

It is not hard to show that t is a proper injection and hence defines a 
homeomorphism (with respect to the induced topology from C^' ) onto 
its image. We conclude with the following global geometric result: 
the smallest affine embedding Xr — ^ is when N = !l!_t^^ which 
follows from the fact that {ti, ...,tN^} is a mmzma/ generating set. 

In Section [5] we discussed that X3 is a branched double cover of C^. 
In general, 

Xr = Spec^„^. (C[ti, ...,tN^]/3r) 

where 3r is an ideal generated by | (^(3)^ + (3) j + ^^(4) polynomials 
(see |Dre03j l 
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